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ABSTRACT 

Using all-sky maps obtained with CORE/DIRBE, we reanalyzed the diffuse sky brightness 
at 1.25 and 2.2/im, which consists of zodiacal light, diffuse Galactic light (DGL), integrated 
starlight (ISL), and isotropic emission including the extragalactic background light. Our new 
analysis including an improved estimate of the DGL and the ISL with the 2MASS data showed 
that deviations of the isotropic emission from isotropy were less than 10% in the entire sky at high 
Galactic latitude (|6| > 35°). The result of our analysis revealed a significantly large isotropic 
component at 1.25 and 2.2/im with intensities of 60.15 ± 16.14 and 27.68 ± 6.21 nWm _2 sr _1 , 
respectively. This intensity is larger than the integrated galaxy light, upper limits from y-ray 
observation, and potential contribution from exotic sources (i.e., Population III stars, intrahalo 
light, direct collapse black holes, and dark stars). We therefore conclude that the excess light 
may originate from the local universe; the Milky Way and/or the solar system. 

Subject headings: cosmic background radiation — dust, extinction — infrared: ISM — infrared: stars 
— scattering — zodiacal dust 


1. INTRODUCTION 

1.1. Extragalactic Background Light in 
the Near-Infrared 

Extragalactic Background Light (EBL) in the 
near-infrared (IR) supposedly comprises inte¬ 
grated light emitted from galaxies, quasars, and 
possible particle decay. Hence, the near-IR EBL 
is a potentially important physical indicator of 
star formation history and the unknown radiation 
processes throughout the history of the universe. 


The lower limit of the near-IR EBL is the 
brightness of the integrated galaxy light (IGL) 
derived from the deep galaxy counts, such as 
those detected in the Hubble Deep Field (HDF) 
by Madau & Pozzetti (2000) or the Subaru Deep 
Field (SDF) by Totani et al. (2001). On the other 
hand, the upper limit of the EBL is estimated from 
observations of high-energy y-ray sources, assum¬ 
ing the property of intrinsic spectra of the objects 
(e.g., Aharonian et al. 2006, Albert et al. 2008, 
Meyer et al. 2012). These y-rays interact with the 
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EBL photons to create positron-electron pairs. At 
1-2 pm, the results of these two methods converge 
to the same intensity (i.e., ~ 10-20 nWm _2 sr _1 ). 

Direct measurement of the EBL is hampered 
by the intense foreground emission, contributed 
by airglow, zodiacal light (ZL), and integrated 
starlight (ISL). In previous studies, the ZL and 
the ISL were subtracted from the sky brightness 
measured by the space telescope. According to 
Matsumoto et al. (2005), who reported on In¬ 
frared Telescope in Space ( IRTS ), and Wright & 
Reese (2000), Cambresy et al. (2001), and Lev- 
enson et al. (2007), who analyzed the Diffuse In¬ 
frared Background Experiment (DIRBE) aboard 
the Cosmic Background Explorer ( COBE ) satel¬ 
lite, the intensity of the EBL in the near-IR can 
be several times higher than that of the IGL and 
that estimated from the high-energy y-ray obser¬ 
vations. If these findings are true, we require addi¬ 
tional sources other than normal galaxies. The ex¬ 
cess light might originate from distant Population- 
Ill (Pop-Ill) stars which cannot be spatially re¬ 
solved by recent observations. However, Inoue et 
al. (2013) calculated the theoretical contribution 
of light from Pop-Ill stars to the EBL, and showed 
that it is smaller than the IGL contribution by 2-3 
orders of magnitude. Therefore, if the excess emis¬ 
sion is real, we must seek other candidate sources. 

1.2. Diffuse Galactic Light 

Previous studies have revealed that diffuse 
Galactic light (DGL) comprises starlight scattered 
off by the interstellar dust grains (e.g., Elvey & 
Roach 1937, Henyey & Greenstein 1941, van de 
Hulst & de Jong 1969, Mattila 1979). The DGL 
contains information on the size distribution of 
interstellar dust grains and the interstellar radi¬ 
ation field (ISRF) that illuminates them. In the 
optically thin limit, the intensity of the far-IR 
100 pm emission is expected to be proportional 
to the DGL intensity (Brandt & Draine 2012). 
Therefore, the DGL has been quantitatively ana¬ 
lyzed by correlating the diffuse optical light with 
the far-IR emission (e.g., Laureijs et al. 1987, 
Guhathakurta & Tyson 1989, Paley et al. 1991, 
Zagury et al. 1999, Matsuoka et al. 2011, Brandt 
& Draine 2012, Ienaka et al. 2013). Although 
the DGL is worthy of study, it constitutes a fore¬ 
ground emission in the EBL measurements, and 
must therefore be removed before analyzing the 


EBL. 

Thus far, the DGL and its contribution to the 
total sky brightness have not been quantified in 
the near-IR. Leinert et al. (1998) suggested that 
the near-IR DGL is limited to low Galactic lati¬ 
tudes (|&| < 5°), where the dust column is suffi¬ 
ciently dense to enhance the intensity of the scat¬ 
tered light. In contrast, Arai et al. (2015) derived 
the mean DGL spectrum at 0.95-1.65 pm using 
the low-resolution spectrometer (LRS) on the Cos¬ 
mic Infrared Background ExpeRiment (CIBER) 
in several local regions of high Galactic latitude 
(|6| > 30°). Although their result is consis¬ 
tent with the DGL spectra modeled by Brandt & 
Draine (2012), whether their result is applicable 
to the general wide field of the sky is not clari¬ 
fied. In addition, Tsumura et al. (2013b) derived 
the DGL spectrum at relatively low Galactic lat¬ 
itudes (5° < |6| < 15°) at 1.8-5.3 pm from data 
collected in the low-resolution prism spectroscopy 
mode of the Infra-Red Camara (IRC) onboard the 
AKARI satellite. To elucidate the general proper¬ 
ties of the DGL and the isotropy of the EBL, these 
studies must be supplemented by measurements of 
the near-IR DGL and its contribution to the sky 
brightness over a wide field. 

1.3. Purpose of the Present Work 

To measure the EBL and DGL at 1 -2 pm, we 
analyze data acquired by the DIRBE aboard the 
COBE spacecraft. For previous studies of the dif¬ 
fuse IR components, the DIRBE observed the all 
sky from 1.25 to 240 pm in 10 bands. These dif¬ 
fuse components included the ZL, thermal emis¬ 
sion from interstellar dust, and the EBL (Hauser 
et al. 1998, Kelsall et al. 1998, Arendt et al. 1998, 
Dwek et al. 1998). Arendt et al. (1998) sub¬ 
tracted the contribution of the ISL from the sky 
brightness using the DIRBE Faint Source Model 
(FSM), which is based on the Wainscoat et al. 
(1992) and Cohen (1993, 1994, 1995) “SKY” mod¬ 
els. As shown in Table 4 of Arendt et al. (1998), 
wherein no entries appear at 1.25 and 2.2 pm, the 
intensity of the 100 pm emission is not correlated 
with diffuse light at those wavelengths. However, 
according to Leinert et al. (1998), the Galactic 
component observed by the DIRBE in these bands 
undoubtedly contains a scattered light contribu¬ 
tion. The missing DGL may have been caused by 
the poor precision of the DIRBE FSM which does 
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not reproduce the actual astrometry and photom¬ 
etry of Galactic stars. 

Since its release, the Two Micron All-Sky Sur¬ 
vey (2MASS) Point Source Catalog (PSC) has 
been used for the starlight subtraction in sev¬ 
eral measurement studies of the near-IR EBL 
(e.g., Wright & Reese 2000, Gorjian et al. 2000, 
Cambresy et al. 2001, Levenson et al. 2007). 
However, because their analyses were limited to 
small regions of low far-IR intensity at high Galac¬ 
tic latitude (|6| > 40°), these authors ignored the 
DGL contribution. 

The present paper reanalyzes the all-sky map 
created by DIRBE for the purpose of evaluating 
the DGL at 1.25 and 2.2 /im and measuring the 
EBL. Calculating the contribution of the ISL col¬ 
lected by the 2MASS PSC over a wide field of 
high Galactic latitudes (|6| > 35°), which includes 
both low and high 100 fim. emission intensity fields, 
we find a positive linear correlation between the 
100 /im emission and the diffuse near-IR light at 
both 1.25 and 2.2/mr. This means that the near- 
IR DGL which was ignored in previous DIRBE 
analyses, is extracted even at high Galactic lati¬ 
tudes. In fields with small DGL components, sub¬ 
tracting the DGL from the isotropic emission does 
not remove the excess brightness against the IGL 
at 1.25 and 2.2 /an, consistent with the previous 
studies. 

The remainder of this paper is organized as 
follows. Section 2 briefly describes the analyzed 
DIRBE data, and Section 3 shows how we esti¬ 
mate the contribution of each near-IR component. 
In this section, the total sky brightness is decom¬ 
posed into the ZL, DGL, ISL, and isotropic emis¬ 
sion components by a \ 2 minimum analysis. Sec¬ 
tion 4 presents the fitting and evaluates the un¬ 
certainty in each component. Section 5, compares 
our fitting results with those of other studies. A 
summary is presented in Section 6. 

Throughout this paper, the surface brightness 
is expressed in nWm' 2 sr _1 or MJy sr _1 . The con¬ 
version formula between these units is 

vl v (nWm _2 sr _1 ) = [3000/A (/zm)] I v (MJysr -1 ). 

( 1 ) 


2. DATA; DIRBE 

DIRBE was primarily designed to search for 
the isotropic IR EBL and to measure its en¬ 
ergy distribution. The cryogenic operation of 
DIRBE was implemented from 1989 November 24 
to 1990 September 21. During these 10 months, 
the sky was observed in 10 bands, from 1.25 /jm 
to 240 /an. The DIRBE instrument was designed 
to make accurate absolute sky-brightness measure¬ 
ments, with a stray light rejection of less than 
1 nWm _2 sr _1 (Magner 1987) and an absolute gain 
calibration uncertainty of 3.1% at 1.25 and 2.2 /nn 
(Hauser et al. 1998). Consequently, the all-sky 
maps at IR wavelengths were created with ~ 0.7° 
beam size. 

Since part of the present study evaluates the 
scaling factor of the DIRBE ZL model (Kelsall 
et al. 1998, hereafter called the “Kelsall model”) 
against the DIRBE data themselves, as described 
in subsection 3.1, we use solar elongation (e) = 90° 
maps from which the ZL is not subtracted. At 
each pixel, the e = 90° maps provide both the 
sky coordinates and the observation date which 
are needed to run the Kelsall model. In contrast, 
Zodi-Subtracted Mission Average (ZSMA) maps 
used in the previous studies (Arendt et al. 1998, 
Cambresy et al. 2001) provide only the sky coordi¬ 
nates at each pixel. Therefore, the Kelsall model 
cannot be used any more in the analysis of the 
ZSMA maps. For this reason, we use the e = 90° 
maps in the present analysis. 

In principle, DIRBE viewed every celestial line 
of sight through the zodiacal cloud at 90° so¬ 
lar elongation once every 6 months; that is, once 
or twice during the 10-month cryogenic mission. 
From the e = 90° maps, we can obtain the IR 
intensity of each wavelength at each pixel by in¬ 
terpolating the observations made at various times 
at e close to 90°. In the following analysis, we use 
the e = 90° maps created through 6 months of 
observation, starting from 1989 January 1. These 
maps cover almost all of the sky. 

Panels (a) and (a’) of Figure 1 illustrate the 
e = 90° map at 1.25 and 2.2 /im, respectively, 
on the Mollweide projection that is reprojected 
from the original “ COBE Quadrilateralized Spher¬ 
ical Cube” (CSC) projection adopted in DIRBE 
products. The CSC projection is an approxi¬ 
mately equal-area projection that projects the ce- 
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et al. 2013, Arai et al. 2015). In the optically 
thin region, the extinction of the DGL is small 
and the correlation is reportedly linear, consis¬ 
tent with theoretical expectation (see Brandt & 
Draine 2012). Theoretically, a linear correlation 
is expected because optically thin fields dominate 
at high Galactic latitudes in the 1.25 and 2.2 /im 
bands. 

Here, we adopt the diffuse 100 pm emission map 
created by Schlegel et al. (1998), hereafter “SFD”, 
which is widely used in the correlation analyses. 
To match the spatial resolutions of the SFD and 
DIRBE maps, we apply an 8 x 8 pixel binning to 
the SFD map. The DGL term Fi(DGL) is then 
defined as 

Fi(DGL) = biF W0 (4) 

where bi is a free parameter and Fjoo is the inter¬ 
stellar 100 /im intensity, defined as follows; 


lestial sphere onto an inscribed cube. In the 
DIRBE convention, each cube face is divided into 
256 x 256 pixels; thus, all-sky maps have 256 2 x 
6 = 393216 pixels. The side of each pixel is ap¬ 
proximately 0.32°. The following analysis is per¬ 
formed on the CSC projection maps. In this pa¬ 
per, we use the e = 90° maps and the beam 
profile maps, available at the DIRBE website 
“lambda.gsfc.nasa.gov/product / cobe/”. 

3. ANALYSIS 

3.1. Model of the Sky Brightness 

The intensity Fj(Obs) of the DIRBE e = 
90° map is evaluated by the brightness model 
Fi (Model), where the subscript refers to one 
of two bands (1.25 or 2.2 pm). Within these bands, 
the sky brightness is assumed as a linear combi¬ 
nation of four components: the ZL, DGL, ISL, 
and isotropic emission including the EBL. The 
Fi (Model) is thus given by 

F t (Model) = Fi(ZL) + F,(DGL) + F;(ISL) + Fff Iso) 

( 2 ) 

where Fj(ZL), Fi(DGL), Fi(ISL), and Fi(Iso) de¬ 
note the intensities of the ZL, DGL, ISL, and 
isotropic emission, respectively. These four terms 
are modeled as follows. 

3.1.1. Zodiacal Light 

The ZL term F./ZL) is defined as 

Fi(ZL) = OiFj(Kel) (3) 

where ai is a free parameter and Fi (Kel) is the ZL 
brightness estimated by the Kelsall model (Kelsall 
et al. 1998). The Kelsall model is a parameterized 
physical model fitted to the time variation of the 
sky brightness measured by the DIRBE. To eval¬ 
uate the scaling factor of the Kelsall model versus 
the DIRBE data, we adopt the free parameter ai. 

If the Kelsall model completely reproduces the sea¬ 
sonal variation of the ZL brightness observed by 
the DIRBE, the parameter a* will equal 1.0. 

3.1.2. Diffuse Galactic Light 

In previous studies of DGL measurements in 
the optical and near-IR, the intensity of the 
100 pm emission was correlated with that of the 
diffuse light (e.g., Matsuoka et al. 2011, Ienaka 


„ f Fsfd — 0.8MJysr 1 (Fsfd > 0.8MJysr *) 
= { O OMJysr-' (Fsfd 5 0.8MJysf) 

(5) 

In the above expression, Fsfd is the 100 pm in¬ 
tensity of the SFD map. Lagache et al. (2000) re¬ 
ported the EBL at 100 pm as 0.78 ± 0.21 MJy sr _1 . 
Matsuoka et al. (2011) correlated the intensi¬ 
ties of the SFD map and optical diffuse light ob¬ 
served by Pioneer 10/11, and revealed a clear 
break around 0.8 MJy sr -1 . Based on these re¬ 
sults, we assumes 0.8MJysr _1 for the EBL at 
100 pm, and subtract this amount from the region 
of Fsfd > 0.8MJysr _1 on the SFD map to obtain 
the 100 pm brightness associated with the inter¬ 
stellar dust. In the region of Fsfd < 0.8MJysr _1 , 
we set Fjoo = O.OMJysr -1 . 


3.1.3. Integrated Starlight 

To estimate the ISL of each region, we created 
integrated brightness maps of the 2MASS sources. 
The 2MASS PSC contains the photometry of ap¬ 
proximately 470,000,000 objects covering 99.998% 
of the sky, with accurate detections below the com¬ 
pleteness limits J = 15.8 and K s = 14.3 mag 
(Skrutskie et al. 2006, Cutri et al. 2006)0 

To convert the magnitudes of the 2MASS 
sources into DIRBE flux densities, we require the 
zero magnitude. As discussed by Cambresy et al. 


1 The 2MASS PSC are available at the website 
‘ http: //www.ipac.caltech.edu/2mass/releases/allsky/doc/explsup. html ’ 
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(a) 1.25um DIRBE (a ) 2.2um DIRBE 



(b) 1.25um 2MASS (|b| > 20 deg) (b’) 2.2um 2MASS (|b| > 20 deg) 



Fig. 1. — Mollweide projections of the all-sky maps at 1.25 (left panels) and 2.2 fim (right panels) in Galactic 
coordinates with the Galactic center in the middle. Panels (a) and (a’) illustrate the full sky DIRBE e = 90° 
intensity maps. The “S” shape depicts the ecliptic plane. Panels (b) and (b’) are integrated brightness maps 
of the 2MASS point sources at high Galactic latitudes (|&| > 20°), created as described in subsection 3.1.3. 
In the present analysis, parts of these maps (|&| > 35°) are used. Panels (c) and (c’) show the intensity 
maps of the DIRBE FSM used by the DIRBE team (Arendt et al. 1998), for comparison with the integrated 
2MASS sources (i.e., panels (b) and (b’)). Each map is arbitrarily scaled. 
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(2001), the difference between the filters of the 
DIRBE and the 2MASS is negligible compared 
with the other uncertainties; hence, they need not 
be corrected. In addition, Levenson et al. (2007) 
correlated the integrated brightness of the 2MASS 
PSC sources against the DIRBE intensity in 40 
high Galactic latitude regions. They reported 
common zero magnitudes at 1.25 and 2.2 /im of 
1467 and 540 Jy, respectively. Accordingly, these 
values are adopted as the zero magnitudes in the 
following analysis. 

We must also apply the DIRBE beam at 1.25 
and 2.2 fin r to the 2MASS sources. An effec¬ 
tive DIRBE beam for acquiring a daily map at 
1.25 fim is illustrated in panel (a) of Figure 2. This 
beam profile measures the relative response of the 
DIRBE to a point source, and includes the sky 
scanning and data sampling effects. In the present 
analysis, the averaged beam should reflect the ob¬ 
servation period of 6 months rather than the daily 
beam, because the intensity of the e = 90° maps 
is the average of dozens of observations. As illus¬ 
trated in panel (b) of Figure 2, the averaged beam 
shapes for the e = 90° maps are estimated by aver¬ 
aging the daily beam profiles. Similar average pro¬ 
files are obtained in both bands (1.25 and 2.2 /im), 
with full width at half-maximums (FWHMs) of 
~ 1 °. 

Since the beam shapes should not largely de¬ 
pend on the location in the CSC projection sky 
map (COBE DIRBE Explanatory Supplement 
1998), we assume that each 2MASS source isotrop¬ 
ically transfers its flux to the nearest 13 pixels on 
the map according to the averaged beam shape 
(panel (b) of Figure 2). Applying this scheme to 
the 2MASS sources at brightnesses below the com¬ 
pleteness limit (J < 15.8 and K s < 14.3mag), we 
calculate the integrated brightness of each pixel 
at high Galactic latitudes (|6| > 20°). Accord¬ 
ing to the Explanatory Supplement to the 2MASS 
All Sky data Release and Extended Mission Prod¬ 
ucts (Cutri et al. 2006), Galactic extinction below 
Galactic latitudes of 35° renders the stellar color 
redder than the intrinsic color. Since this effect 
attenuate the EBL and disrupt the linear com¬ 
bination of the fitting model (Equation (2)), we 
limit the following analysis to the high Galactic 
latitude regions (|6| > 35°). Panels (b) and (b’) 
of Figure 1 are integrated brightness maps of the 
2MASS sources at |6| > 20° at 1.25 and 2.2 /im, 


respectively. For comparison, the FSM map used 
by the DIRBE team (Arendt et al. 1998) is also 
shown (see panels (c) and (c’) of Figure 1). In 
contrast to the FSM maps, wherein the surface 
brightness smoothly changes across the sky, the 
2MASS-derived maps show clear fluctuations re¬ 
flecting the astrometry and photometry of the 
actual sources. 

In terms of the integrated brightness of the 
2MASS sources (Ei(2MASS); see Figure 2), the 
total ISL term E,(ISL) is defined as 

Fi(ISL) = CjEj(2MASS) (6) 

where Ci is a free parameter representing the inte¬ 
grated brightness of stars fainter than the limiting 
magnitude of the 2MASS. This formula assumes 
that the integrated brightness of the fainter stars 
and the brighter sources below the 2MASS detec¬ 
tion limit have the same spacial distribution. In 
previous studies using the 2MASS data for star 
subtraction (e.g., Cambresy et al. 2001, Wright 
2001), the analyzed region was sufficiently small 
to assume isotropic ISL of fainter stars; thus the 
contributions of fainter stars were subtracted by 
star-counts models (e.g., Jarrett, SKY model). In 
contrast, the present analysis covers a wide field 
of the high-latitude Galactic sky, where the ISL of 
fainter stars and brighter sources should have the 
same spatial distribution. This justifies our use of 
Equation (6), which is free from the uncertainties 
introduced by the star-counts model. The system¬ 
atic features of this simple model are discussed in 
subsection 5.2. 

The 2MASS PSC should contain faint galax¬ 
ies that are not resolved as extended sources. 
The contributions of these faint objects should 
be isotropic in the sky and should be included in 
the EBL. Wright (2001) estimated that galaxies 
with magnitudes K s < 14.3 mag contribute ap¬ 
proximately order of 0.12 and 0.14nWm~ 2 sr _1 at 
1.25 and 2.2 /an, respectively. To derive the EBL 
intensity, we apply these small corrections to the 
isotropic term di after the fitting procedure. 

3.1.4- Isotropic Emission 

Since the isotropic emission is assumed to be 
independent of the region, the term Fi(Iso) is de¬ 
fined as 

^(Iso) = d t (7) 
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Fig. 2.- Profiles of the DIRBE beams. Panel (a) is the DIRBE beam for daily maps (FWHM ~ 0.7°) 
at 1.25/zm, available as “Beam Profile Maps”, in the COBE/DIRBE website. Panel (b) is an averaged 
beam (FWHM ~ 1°), created by averaging the effective beams of the daily maps shown in panel (a). The 
integrated brightness maps in panels (b) and (b’) of Figure 1 were created by applying this averaged beam 
to each 2MASS source. 


where di is a free parameter. 

3.2. Fitting 


To determine the parameters bi, Ci, and 
di, we minimize the following y 2 function in each 
band; 


Now the model brightness Fj (Model) for Fj(Obs) 
is given by 

F % (Model) = Fi( ZL) + F,(DGL) + F,(ISL) + F t ( Iso) 

( 8 ) 

= Oi Fi (Kel) + b x F 100 + c, Fi (2MASS) + d t (9) 

Prior to fitting, we should remove pixels that 
might perturb the analysis. To suppress the large 
photometric uncertainty of bright stars, we mask 
the pixels around stars brighter than J = 5 and 
K s = 4 mag on the CSC projection map at 1.25 
and 2.2 fix n, respectively. In addition, we blank out 
the circular regions around the Magellanic Clouds 
and probable Galactic extended sources listed in 
the Explanatory Supplement to the 2MASS All 
Sky Data Release and Extended Mission Products 
(Cutrietal. 2006). Furthermore, we select regions 
with F S fd < lOMJysr -1 (where the Galactic ex¬ 
tinction is assumed negligible) and exclude outliers 
by applying 2 sigma clipping to the integrated in¬ 
tensities of the 2MASS sources Fi(2MASS). Ap¬ 
proximately 65% of the total pixels in the 16| > 35° 
region in both bands survive these masking proce¬ 
dures. 


: ^2 [ f *( Qbs ) - (Model)] 2 ^ 

j <J ' 1 


[Fj(Obs) 

j 


a* Ft (Kel) - WF, 00 - c t F t ( 2MASS) - d t } 2 


( 11 ) 

where “j” refers to the pixels used in the fitting. 
The total uncertainty ay in each pixel is calculated 
as follows: 


a- 2 = o-j (Obs) 2 + b 2 af 00 + cj<7i (2MASS) 2 (12) 

where (Obs), ctioo, and o-i(2MASS) are the stan¬ 
dard deviations of the intensities in the e = 90° 
map, the intensities of the 100 fim. emission, and 
the integrated intensity of the 2MASS sources, re¬ 
spectively. We adopt aioo = 0.35MJysr _1 de¬ 
rived by Ienaka et al. (2013). The oy(2MASS) at 
each pixel is calculated identically to the bright¬ 
ness of the 2MASS sources (see subsection 3.1.3): 

oy(2MASS) 2 = [-0.4 (loglO) 10-°- 4 "V TOi F l0 ] 2 

(13) 
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where m*, <J mi , and Fw denotes the magnitude of 
each 2MASS source, the uncertainty of this magni¬ 
tude, and the zero magnitude derived by Levenson 
et al. (2007), respectively. Sources lacking a pho¬ 
tometric uncertainty entry in the 2MASS PSC are 
assigned an uncertainty of 0.5 mag. 

4. RESULTS 

4.1. Results of the Decomposition 

The parameters determined by the fitting at 
1.25 and 2.2 pm in the |6| > 35° region are sum¬ 
marized in Table 1. Owing to the large sample 
size (over 100,000 points), the statistical uncer¬ 
tainty in each parameter is very much smaller 
than the determined value. In Figure 3, the 
sky brightness obtained by the DIRBE observa¬ 
tions is decomposed into the ZL, DGL, ISL, and 
isotropic emission according to the linear combina¬ 
tion model (Equation (2)). Filled circles represent 
the weighted means of the points within an arbi¬ 
trary x-direction bin. In further discussion, these 
weighted means will be assumed as representative 
values. 

As shown in panels (b) and (b’) of Figure 3, 
the diffuse near-IR light is positively linearly cor¬ 
related with the interstellar 100 /im emission at 
both 1.25 and 2.2 /im, and the correlations are 
significant. This indicates that the DGL compo¬ 
nent certainly exists, even at high Galactic lati¬ 
tudes (|6| > 35°). Moreover, the linear correlation 
continues through the low to high 100 /jm inten¬ 
sity region (F 10 o < 9MJysr _1 ), indicating that 
the 100 pm emission also well-tracks the DGL at 
these wavelengths. This trend was discovered ow¬ 
ing to the wide sky coverage of the DIRBE maps, 
which have wide dynamic range of the 100 pm in¬ 
tensity. In conclusion, by combining the precise 
star subtraction from the 2MASS PSC data with 
wide-held coverage of the sky, we can find the near- 
IR DGL at 1.25 and 2.2 pm even at low column 
density. 

As shown in panels (a), (a’), (c), and (c’) of Fig¬ 
ure 3, the ZL and ISL are also decomposed from 
the sky brightness with high linearity. Therefore, 
the assumed linear combination model of the sky 
brightness is appropriate for our purpose. 

Although the residuals F)(Obs) — F t (Model) 
in each panel of Figure 3 appears to be func¬ 
tions of Fj(Kel), Fioo, and Fj(ISL), they deviate 


from the best-fit line by no more than ±10 and 
±5nWm -2 sr -1 at 1.25 and 2.2 pm, respectively. 
As shown in Table 2, these deviations are within 
±2% of the typical DIRBE intensity Fi(Obs) in 
both bands. Therefore, the large sample size of 
high quality DIRBE data has enabled a very pre¬ 
cise analysis. The possible origins of the sys¬ 
tematic features in the residuals are discussed in 
section 5.2. 

4.2. Uncertainty Estimation of the Param¬ 
eters 

The statistical uncertainty of the parameters 
determined in the minimum y 2 analysis is not the 
only uncertainty in each component. Other un¬ 
certainties include the parameter variation among 
different regions, the absolute gain of the DIRBE, 
the uncertainty in the faint galaxies in the 2MASS 
PSC, and the systematic uncertainty in the Kelsall 
model. 

4-2.1. Regional Parameter Variations 

To understand the parameter variation among 
different regions with similar dynamic ranges of 
each component, we divide |6| > 35° region into 
6 Galactic longitude fields, i.e., 0° < l < 60°, 
60° < l < 120°, 120° < l < 180°, 180° < l < 240°, 
240° < l < 300°, and 300° < l < 360°. In each 
region, we apply the fitting procedure described 
in subsection 3.2. The results for each region are 
shown in Table 1. Because each region contains 
over 10,000 points, the statistical uncertainty in 
this analysis remains small. Figure 4 presents 
the parameters obtained in the 6 regions and in 
the \b\ > 35° region as functions of Galactic lon¬ 
gitude. Reasonably, the parameters determined 
by the fitting in the 6 regions are randomly dis¬ 
tributed around the parameters determined in the 
161 > 35° field, showing some degree of scatter. 
The standard deviation of the parameter values in 
the 6 regions is adopted as a conservative uncer¬ 
tainty and is listed for each parameter in the row 
“Scatter” in Table 3. 

4-2.2. Uncertainty in the Absolute Gain of 
DIRBE 

Hauser et al. (1998) reported an uncertainty 
of 3.1% in the absolute gain of the DIRBE at 
1.25 and 2.2 pm. The Kelsall model was devel- 
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Fig. 3.— Fitting results at 1.25 and 2.2/rm (|6| > 35°). Top panels (a) and (a’) plot 

Fi(Obs) — Fi(DGL) — Fi(ISL) (i.e., a,iFi( Kel) + di ) versus Fi(Kel); center panels (b) and (b’) plot 

Fi(Obs) — Fi( ZL) — Fi(ISL) (i.e., biFioo + di) versus Fioo, and bottom panels (c) and (c’) plot 

Fi(Obs) — Fi( ZL) — Fi(DGL) (i.e., CiFi(2MASS) + di) versus Fi(2MASS). The red lines are the results of 

the best-fit parameters. The middle and bottom parts of each panel plot the residuals Fi( Obs) — Fi (Model) 
and the number of pixels, respectively, as functions of Fi(Kel) (top), Fioo (center), and Fi(2MASS) (bot¬ 
tom). The filled circles and error bars represent the weighted means and the weighted standard errors of the 
sample within the arbitrary x-direction bin. 
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Table 1: Results of the fitting in each region 


Band 

Region 

Number of 

CLi 

vbi 

Ci 

vdi 

(H 

(deg) 

pixels 

(dimensionless) 

(nWm -2 MJy -1 ) 

(dimensionless) 

(nWm -2 sr -1 ) 

1.25 

|6| > 35 

116578 

1.0087 ±0.0001 

4.79 ± 0.02 

1.0238 ±0.0003 

60.03 ±0.08 

2.2 

|6| > 35 

119394 

1.0450 ± 0.0002 

1.49 ±0.01 

1.0333 ±0.0004 

27.54 ±0.04 

1.25 

|6| > 35, 0 < l < 60 

18557 

1.0017 ±0.0004 

8.67 ±0.06 

1.0206 ±0.0009 

64.51 ±0.25 

2.2 

6 > 35, 0 < l < 60 

18852 

1.0350 ± 0.0006 

2.79 ± 0.03 

1.0389 ±0.0011 

27.90 ±0.13 

1.25 

|6| > 35, 60 < l < 120 

19714 

0.9956 ± 0.0003 

1.91 ±0.05 

1.0276 ±0.0008 

64.25 ± 0.18 

2.2 

6 > 35, 60 < l < 120 

20309 

1.0286 ± 0.0004 

0.94 ± 0.03 

1.0409 ±0.0010 

29.15 ±0.10 

1.25 

|6| > 35, 120 <1 < 180 

19582 

0.9993 ± 0.0003 

5.60 ±0.04 

1.0236 ±0.0009 

68.51 ±0.20 

2.2 

\b\ > 35, 120 <1 < 180 

20213 

1.0534 ±0.0005 

1.36 ±0.02 

1.0288 ±0.0011 

27.79 ±0.10 

1.25 

6 > 35, 180 <1 < 240 

20275 

1.0229 ±0.0003 

4.07 ±0.04 

0.9986 ± 0.0008 

55.97 ±0.20 

2.2 

\b\ > 35, 180 <1 < 240 

20748 

1.0492 ± 0.0005 

1.26 ±0.02 

1.0196 ±0.0011 

28.88 ±0.11 

1.25 

\b\ > 35, 240 <1 < 300 

19516 

1.0224 ±0.0002 

0.64 ± 0.06 

1.0223 ±0.0008 

53.92 ±0.19 

2.2 

\b\ > 35, 240 < l < 300 

20011 

1.0593 ± 0.0004 

-0.07 ±0.04 

1.0306 ±0.0010 

26.12 ±0.10 

1.25 

6 > 35, 300 <1 < 360 

18934 

1.0030 ± 0.0003 

5.58 ±0.05 

1.0421 ±0.0008 

58.62 ±0.23 

2.2 

6 > 35, 300 <1 < 360 

19261 

1.0372 ± 0.0005 

2.04 ±0.03 

1.0524 ±0.0010 

25.87 ±0.12 


Note. - The symbols in the column headings are defined in Section 3. 

Error in each component is the statistical uncertainty derived by the fitting. 


Table 2: Typical intensity of each component determined by the fitting at |6| > 35° 


Component (nWm "^sr i ) 

1.25 pm 

2.2 pm 

Fi{ ZL) = ai F t { Kel) 

544 ± 199 

208 ± 74 

F t (DGL) = biF W0 

4.87 ±6.06 

1.50 ± 1.87 

Fi( ISL) = c I E l (2MASS) 

175 ± 80 

66.2 ±34.7 

Fi (Iso) = d t 

60.03 ±0.08 

27.54 ±0.04 

Fi( Obs) 

787 ± 220 

304 ± 87 


Note. - Except for F*(Iso), each component is represented 
by its average and standard deviation of the samples in the |6| > 35° region, 
where the regions around the 2MASS sources brighter than J = 5 and K a = 4 mag 
at 1.25 and 2.2 pm are masked, respectively. 
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Fig. 4.— Parameter variation among the 6 sampled regions as functions of the Galactic longitude. Black 
and red symbols are the parameter values determined at 1.25 p,m and 2.2/im, respectively. Panels (a), (b), 
(c), and (d) plot the parameters at, bi, Ci and di, respectively. Circles represent the results in each of the 
6 regions and at |6| > 35°. Horizontal error bars indicate the range of the region. Horizontal dashed lines 
represent the average values in the 6 regions. 


oped to match the photometric scale of DIRBE 
(Kelsall et al. 1998). The SFD map and the ISL 
of the 2MASS sources are also scaled to the pho¬ 
tometric scale of the DIRBE by the fitting pro¬ 
cess. Then the parameters a*, bi, and c, are un¬ 
affected by the uncertainty in the absolute gain. 
However, this uncertainty influences the parame¬ 
ter di . The value of these uncertainties (assuming 
a percentage contribution of 3.1%) appear in the 
row “Gain” in Table 3. 

4-2.3. Uncertainty Associated with the Faint 
Galaxies in the 2MASS PSC 

As explained in subsection 3.1.3, the 2MASS 
PSC may contain unresolved faint galaxies in ad¬ 
dition to the Galactic stars. Wright (2001) es¬ 
timated that galaxies with K s < 14.3 mag con¬ 
tribute around 0.12 and 0.14nWm _2 sr _1 to the 
isotropic emission at 1.25 and 2.2 /im, respec¬ 
tively. Therefore, we add these corrections to the 
isotropic term d t after decomposing the integrated 
brightness of stars. The contribution of unresolved 
galaxies is also added to the uncertainty of d t and 
is listed in the “Galaxies” row in Table 3. Note 
that these contributions are relatively small. 


4.2.4. Systematic Uncertainty in the Kelsall 
Model 

As reported in Kelsall et al. (1998), the uncer¬ 
tainty in the ZL model is 15 and 6nWm _2 sr _1 at 
1.25 and 2.2 yrm, respectively. These values were 
estimated as the difference between two ZL mod¬ 
els in the north Galactic pole region. These two 
models were equally good in reproducing the ob¬ 
served seasonal variations in the ZL but not in the 
isotropic component. The contributions of these 
uncertainties to the uncertainty in di are listed in 
the row “ZL model” in Table 3. 

4-2.5. Total Uncertainty 

The quadrature sum of the uncertainties in each 
parameter is presented in the row “Quadrature 
sum” in Table 3. In the following Discussion sec¬ 
tion, the parameters of the ZL, DGL, ISL, and 
isotropic emission components are assumed as the 
parameters determined in the |6| > 35° region and 
their errors are assumed as the quadrature sums 
of the uncertainties. 
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Table 3: Uncertainty budget associated with each parameter 


Band (/jm) 

at (dimensionless) 

ubi (nWnr 2 MJy i ) 

d (dimensionless) 

vdt (nWm 2 sr i ) 

1.25 

2.2 

1.25 

2.2 

1.25 

2.2 

1.25 

2.2 

Statistical 

0.0001 

0.0002 

0.02 

0.01 

0.0003 

0.0004 

0.08 

0.04 

Scatter 

0.012 

0.012 

2.88 

0.97 

0.014 

0.011 

5.66 

1.37 

Gain 

— 

— 

— 

— 

— 

— 

1.86 

0.85 

Galaxies 

— 

— 

— 

— 

— 

— 

0.12 

0.14 

ZL model 

— 

— 

— 

— 

— 

— 

15 

6 

Quadrature sum 

0.012 

0.012 

2.88 

0.97 

0.014 

0.011 

16.14 

6.21 

Result ( 6 > 35°) 

1.0087 

1.0450 

4.79 

1.49 

1.0238 

1.0333 

60.03 

27.54 


Note. - Symbols in the column headings are defined in Section 3. 


5. DISCUSSION 

5.1. Interpretation of the Determined Pa¬ 
rameters 

As shown in Table 3, the parameter ai at 
1.25 /im (determined as 1.0) lies within the un¬ 
certainty limits, indicating that the Kelsall model 
well-reproduces the time variation of the sky 
brightness measured by the DIRBE at 1.25 pm. 
The parameter cp at 2.2 pm exceeds 1.0 by ap¬ 
proximately 4%. On average, a 4% variation in the 
Kelsall model corresponds to 8 nWnr _2 sr , which 
is slightly larger than the claimed systematic un¬ 
certainty in the model (6nWm _2 sr _1 ). This sug¬ 
gests that the Kelsall model underestimates the 
ZL intensity in this band. For one thing, in the 
fitting procedure, the Kelsall model was sampled a 
sky pixel every ~ 5° or ~ 10° as a spatial grid (not 
used the all pixels) to avoid the excessive compu¬ 
tational requirements. In addition, comparing the 
parameter values at 1.25 and 2.2 pm, some of the 
parameters determined in the Kelsall model, espe¬ 
cially the phase function parameter Gy 2 and the 
Albedo A 2 , seem anomalous 2.2pm [see Table 2 
of Kelsall et al. (1998)]. Specifically, at 2.2pm, 
Gy 2 is 3 times smaller than at 1.25 pm and A 2 is 
unnaturally larger than that at 1.25 pm. These 
results may change the spatial distribution of the 
ZL brightness in the Kelsall model, and may also 
explain why ai deviate from 1.0 at 2.2 pm. 

The uncertainty in the parameter bi is domi¬ 
nated by scatter among the different regions and 
exceeds 50% of the result. This large error is 
attributed to the low typical brightness of the 


DGL (1-2 orders of magnitude fainter than the 
other components; Table 2). In that situation, if 
the SFD intensity spatially correlates with that 
of the Kelsall model or the ISL to some extent, 
some of the DGL might be absorbed by the other 
components in the fitting process. Investigating 
this possibility is beyond the scope of the present 
study; instead, we conservatively estimate the un¬ 
certainty in the DGL as the scatter in the 6 re¬ 
gions. Remarkably, the present analysis identified 
the DGL despite its much lower brightness than 
that of the other components, by virtue of the well- 
calibrated all-sky maps of the DIRBE. The DGL 
results in the optical and near-IR, determined in 
the present and previous studies, are compared in 
subsection 5.3. 

The parameter Q exceeds 1.0 by 1-4% in both 
bands. Assuming that the ISL of the sources 
brighter and fainter than the 2MASS detection 
limit have the same spatial distribution, this ex¬ 
cess is presumably contributed by the fainter stars. 

This result also indicates that the zero magnitude 
derived by Levenson et al. (2007) is appropri¬ 
ate for converting the 2MASS magnitude to the 
DIRBE flux in the present study. 

The determined parameter di is discussed in 
section 5.4. 

5.2. Dependence of the Residuals on Galac¬ 
tic and Ecliptic Latitude 

Figure 5 illustrates the residuals Obs) — Fi (Model) 
derived from the fitting in the |6| > 35° region as 
functions of Galactic latitude b and ecliptic lati¬ 
tude /3. In general, the dependence of the residuals 
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on Galactic latitude traces the ISL or the DGL, 
whose intensities are also functions of Galactic 
latitude. On the other hands, the dependence on 
ecliptic latitude is expected to measure the accu¬ 
racy of the ZL model. 

As shown in panels (a) and (a’) of Figure 5, 
the residuals tend to increase toward low Galac¬ 
tic latitudes, suggesting that some component is 
missed at these latitudes. Similar trends are ob¬ 
served when the residuals are plotted against the 
integrated intensity of the 2MASS sources (panels 
(c) and (c’) in Figure 3). In this case, the residuals 
increase toward regions of higher intensity, where 
the data of the lower Galactic latitude fields are 
more dominant. 

The phenomenon might stem from the contri¬ 
bution of stars with no entry in the 2MASS PSC, 
possibly because they were masked by their near¬ 
est bright sources. According to the Explana¬ 
tory Supplement to the 2MASS All-Sky Data Re¬ 
lease and Extended Mission Products (Cutri et 
al. 2006), masking around bright stars can filter 
faint sources from the detection process. Although 
the masked area in the all-sky averages to 0.25% 
and 0.43% at 1.25 and 2.2 ^m, respectively, the 
fraction of such regions tends to increase toward 
lower Galactic latitudes as the number density of 
bright sources increases. In addition, the 2MASS 
compensated for saturation caused by bright stars 
by fitting the unsaturated wings of their intensity 
profiles. This suggests that the 2MASS PSC could 
have missed faint stars. 

The simply modeled ISL term Fi(ISL) = aiFi(2MASS) 
might also contribute to the latitude dependence 
of the residuals. If the integrated intensities of the 
bright and faint stars (below the detection limit of 
the 2MASS) have different spatial distributions, 
the model assumption is not strictly valid. Al¬ 
though the dominant cause of the features in the 
residuals cannot be determined, the amplitude of 
the residuals is within ±10 and ±5nWm - 2 sr -1 
at 1.25 and 2.2 /tm, respectively. Within these 
ranges, the di terms are certainly isotropic. The 
origin of the residuals’ dependence on b requires 
searching by an all-sky with superior sensitivity 
and spatial resolution to the 2MASS. 

As illustrated in panels (b) and (b’) of Figure 
5, the dependence of the residuals on ecliptic lat¬ 
itude, especially the turbulence near the ecliptic 
plane, may reflect the incompleteness of the Kel- 


sall model. When the residuals are plotted against 
the Kelsall model Fi(Kel) (panels (a) and (a’) of 
Figure 3), distortion appears in the higher inten¬ 
sity region (comprising fields of lower ecliptic lati¬ 
tudes). Cambresy et al. ( 2001 ), who similarly sub¬ 
tracted the ZL using the Kelsall model, reported 
the same trend. These results highlight the dif¬ 
ficulty in applying the ZL model near the eclip¬ 
tic plane, where the distribution of the zodiacal 
dust (including the dust bands and the circumso¬ 
lar ring) becomes complex. 

The effects of these latitude dependences are 
naturally included in the scatter of the fitting 
results among the different regions (Figure 4). 
Therefore, the fluctuations related to Galactic or 
ecliptic latitude are not added to the uncertainty 
budget. 

5.3. Spectrum of the Diffuse Galactic 
Light 

We now discuss the spectrum of the parame¬ 
ter bi in the optical and near-IR. The four col¬ 
ored curves in Figure 6 are synthetic DGL spec¬ 
tra calculated by Brandt & Draine (2012) based 
on two estimates of the ISRF continuum and two 
dust models; namely, the Zubko et al. (2004) 
and Weingartner & Draine (2001) models, here¬ 
after referred to as ZDA04 and WD01, respec¬ 
tively. The WD01 model assumes that the half¬ 
mass grain radius ao .5 (denoting that 50% of the 
mass is contributed by grains with radii a > ao. 5 ) 
is ~ 0.12 /im for both silicate and carbonaceous 
grains (Draine 2011). In the ZDA04 model, grains 
with radii a > 0.2 /im comprises a small propor¬ 
tion of the mass, and the half-mass radius differs 
between carbonaceous grains (ao .5 ~ 0.06 /.an) and 
silicate grains (ao .5 ~ 0.07 /im). Consequently, the 
WD01 creates a redder scattered spectrum than 
the ZDA04 at far-optical and near-IR wavelengths. 
The local ISRF continua are estimated either from 
Mathis et al. (1983) with de-reddening of the orig¬ 
inal ISRF of the MMP83 (see Brandt & Draine 
2012 for details) or from a synthesis model of solar- 
metallicity star populations (Bruzual & Chariot 
2003). The de-reddening correction and Bruzual - 
Chariot model are hereafter referred to as MMP83 
and BC03, respectively. In the latter, the star 
formation rate is proportional to exp(—f/5 Gyr), 
where t denotes the star formation timescale in 
units of Gyr. 


13 





Fig. 5.- Dependence of the residuals Fi( Obs) — Fi (Model) on Galactic and ecliptic latitudes, determined 
by the fitting in the | b\ > 35° region. The upper part of each panel plots Fi( Obs) — .Fj (Model) at 1.25 and 
2.2 ^m as a function of Galactic latitude b (top panels) and ecliptic latitude f3 (bottom panels). Filled circles 
and errors bars denote the weighted means and the errors of the points within arbitrary x-direction bins, 
respectively. The lower part of each panel is a histogram of the number of pixels at each b or /3. 
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Fig. 6.— Parameter bi as a function of wavelength in different analysis. Plotted are the present results 
(filled red circles), the Pioneer 10/11 data (Matsuoka et al. 2011) (open circles), results of a translucent 
cloud at high Galactic latitude (Ienaka et al. 2013) (open triangles), results derived from GIBER (Arai et al. 
2015) (open squares), and results derived from AKARI (Tsumura et al. 2013b) (asterisks connected by the 
solid line). The colored curves are the synthetic spectra of the ratio of DGL to 100/mr intensity calculated 
by the WD01/BC03 model (green dash-dotted curve), the WD01/MMP83 model (green dashed curve), the 
ZDA04/BC03 model (blue dash-dotted curve), and the ZDA04/MMP83 model (blue dashed curve). To fit 
the observed results in the optical, the WD01 and ZDA04 are scaled by 1.7 and 1.9 times, respectively. 
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In the optical region, Figure 6 plots the results 
collected by Pioneer 10/11 (Matsuoka et al. 2011) 
(open circles) and toward a high Galactic latitude 
translucent cloud MBM32 (Ienaka et al. 2013) 
(triangles). The Pioneer 10/11 results were ob¬ 
tained in the same field of the sky as the present 
analysis (|6| > 35°). In the near-IR, Figure 6 
plots our results (filled red circles), the mean of 
6 small regions from GIBER (Arai et al. 2015) 
(squares), and the AKARI results (Tsunrura et al. 
2013b), collected at relatively low Galactic lati¬ 
tudes (5° < |6| < 15°) (asterisks). 

Around 1.25 pm, our result is marginally con¬ 
sistent with the GIBER results, considering the 
uncertainties in the measurements. This suggests 
that results acquired in local regions represent the 
wider region at high Galactic latitudes. At 2.2 /rnr, 
our result is consistent with the AKARI results, 
allowing for the uncertainties. In contrast to our 
results at high Galactic latitude (|6| > 35°), the 
AKARI results were taken at lower galactic lati¬ 
tude (5° < |6| < 15°). Whether the relationship 
between the parameter b t and Galactic latitude 
results from the wide scatter of bi among the dif¬ 
ferent regions (subsection 5.1) is difficult to deter¬ 
mine. 

The green and blue curves show the spectra of 
the scaled ZDA04 and WD01, respectively. Since 
the original ZDA04 and WD01 models underesti¬ 
mate the observed bi by a factor of 2, these models 
are arbitrarily scaled to the optical results (Mat¬ 
suoka et al. 2011, Ienaka et al. 2013) by factors 
of 1.9 and 1.7, respectively. Ienaka et al. (2013) 
suggested two possible explanations for this dis¬ 
crepancy: deficient UV photons in the ISRF and 
underestimation of the assumed albedo of the dust 
grains in the models. Combining the results in 
the optical, GIBER, and the present study, we 
find that the scaled ZDA04 provides a better fit¬ 
ting spectrum than WD01 at 1 — 2 pm. imply¬ 
ing a bluing of the DGL spectrum in this wave¬ 
length range. In contrast, the spectrum derived 
from AKARI significantly exceeds the ZDA04 and 
WD01 spectra at longer wavelengths (> 3 pm), 
possibly because it includes the thermal emission 
at low Galactic latitudes, whereas the ZDA04 and 
WD01 spectra include only the scattered light 
component. At shorter wavelengths (< 3 pm), 
the observed DGL can be well-fitted to the model 
spectra containing the scattered component alone. 


5.4. AN ISOTROPIC EMISSION COM¬ 
PONENT 

We estimate the EBL intensity from the derived 
parameter di. This estimate adds the contribu¬ 
tion of the faint galaxies appearing in the 2MASS 
PSC back to di. Assuming the faint-galaxy con¬ 
tribution estimated by Wright (2001) and the un¬ 
certainty in di (Table 3), the EBL intensity at 
1.25 and 2.2 /im is estimated as 60.15 ± 16.14 and 
27.68 ± 6.21 nWm _2 sr _1 , respectively. 

5.4-1- Isotropy test 

We now discuss the isotropy of the EBL. Con¬ 
sidering only the scatter of di among the different 
regions (Table 3) and disregarding the other un¬ 
certainties, the deviations from isotropy are less 
than 10% and 5% of the determined di at 1.25 
and 2.2 pm, respectively. These isotropies are con¬ 
sistent with the typical relationships between the 
residuals and the Galactic and ecliptic latitudes 
(Figure 5). This suggests that the EBL is isotropic 
within these limits even in the strong foreground 
emission. The EBL isotropy can be usefully ex¬ 
amined in the present study because of the wider 
observation region than in previous studies. 

5-4-2. Comparison with Other Studies 

Figure 7 compares the resultant EBL with 
those of previous studies. Using the DIRBE 
data, Cambresy et al. (2001) derived the EBL 
at 1.25 and 2.2 pm, by subtracting the Galactic 
stars by 2MASS and the ZL by the Kelsall model. 
These authors targeted regions with low intensity 
of the dust emission (DIRBE 240 pm brightness 
I240 < 3MJysr _1 ). In such regions, the expected 
DGL brightness at 1.25 and 2.2 pm is < 7 and 
< 2nWnr~ 2 sr _1 , respectively, assuming the pa¬ 
rameter bi determined in the present study and 
the conversion factor between the 100 and 240 pm 
intensities (1.297; see Table 4 of Arendt et al. 
(1998)). For this reason, Cambresy et al.’s (2001) 
study ignored the DGL. In addition, the -F 100 his¬ 
tograms (panels (b) and (b’) of Figure 3) show 
that regions of lower 100 pm intensity (and DGL 
brightness) dominate in the sky. Therefore, the 
EBL results derived in this study are reasonably 
consistent with those obtained by Cambresy et 
al. (2001), despite the lack of any quantitative 
DGL evaluation in the latter study. However, 
Cambresy et al. (2001) noticed fluctuations in 
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EBL with ecliptic latitude, which they attributed 
to a small DGL component at 1.25 /im [see Figure 

5 of Cambresy et al. (2001)]. 

Using the FSM (Faint Source Model) for the 
starlight subtraction, Hauser et al. (1998) derived 
the EBL intensity at high Galactic and ecliptic lat¬ 
itudes. Although their results are plotted as 95% 
confidence upper limits in Figure 7, their direct 
values are significantly smaller than ours; 33.0±21 
and 14.9 ± 12nWm _2 sr _1 at 1.25 and 2.2 /an, re¬ 
spectively (see Table 2 of Hauser et al. 1998). 
This discrepancy can be explained by the following 
two things related to the ISL evaluation. At first, 
in converting the magnitudes of the sources into 
DIRBE flux densities, the present study adopts 
the zero magnitude of 1467 and 540 Jy at 1.25 
and 2.2 /tin, respectively, but Hauser et al. (1998) 
used higher one, i.e., 1547 and 612.3 Jy at 1.25 and 
2.2 /im, respectively (COBE DIRBE Explanatory 
Supplement 1998). We used the zero magnitudes 
derived by Levenson et al. (2007), who corre¬ 
lated the intensity of the 2MASS-derived ISL with 
that of the DIRBE and corrected the zero magni¬ 
tude to fit the photometric scale of the 2MASS 
to that of the DIRBE. Therefore, the zero magni¬ 
tudes we adopted are suitable to estimate the ISL 
contribution in the DIRBE data. Next, Wright 

6 Reese (2000) suggested that the Wainscoat et 
al. (1992) star-counts model, which is the basis of 
the FSM, overestimates the counts by ~ 10% in 
the 6 < K < 10 range at high Galactic latitudes, 
compared with the 2MASS. This is within the 10- 
15% uncertainty of the FSM, estimated in Arendt 
et al. (1998). Considering these differences asso¬ 
ciated with the ISL estimation, the ISL intensity 
in Hauser et al. (1998) can be higher than that 
in the present study by ~ 15 and ~ 20% at 1.25 
and 2.2 /im, respectively. These percentages cor¬ 
respond to ~ 26 and ~ 13nWm^ 2 sr _1 at 1.25 
and 2.2 pa n, respectively, assuming the ISL inten¬ 
sity derived in the present study (Table2). This 
overestimation of the ISL in Hauser et al. (1998) 
well explains the resultant EBL differences be¬ 
tween Hauser et al. (1998) and the present study 
at both bands. 

In the EBL measurement, the removal of the 
ZL from the sky brightness is controversial, as 
multiple ZL models are available. For instance, 
Gorjian et al. (2000), Wright (2001) and Leven¬ 
son et al. (2007) used the ZL model based on 


Wright (1998), whereas Cambresy et al. (2001) 
and the present study adopted the Kelsall model. 
As noted by Levenson et al. (2007), the ZL inten¬ 
sity at the ecliptic pole at 1.25 and 2.2 /im is ~ 22 
and ~ 5nWm _2 sr _1 lower in the Kelsall model 
than in Wright’s (1998) model, respectively. Con¬ 
sequently, the difference between the two models 
tends to be larger at 1.25 pm than at 2.2 pm. As 
shown in Figure 7, the resultant EBL obtained 
with the Kelsall model can be a few times lower 
than that obtained with the Wright model espe¬ 
cially at 1.25 pm. At 2.2 pm, the results of both 
models converge within their uncertainties. To 
eliminate the uncertainty introduced by the ZL 
model itself and its variation, we must observe out¬ 
side the ZL cloud. 

Note that the present decomposition analysis 
cannot identify where the isotropic emission in¬ 
cluding the EBL comes from. For example, if the 
isotropic component associated with the ZL ex¬ 
ists, it may contribute to the “EBL” called in this 
paper. Actually, the Kelsall model was developed 
to fit to the seasonal variation of the DIRBE sky 
brightness, ignoring the uniform component if ex¬ 
ists. Hauser et al. (1998) also emphasized that 
the Kelsall model cannot uniquely determine the 
true ZL signal; in particular an arbitrary isotropic 
component could be added to the model. 

5-4-3. Implications of the Present Results 

Consequently, the present EBL result still re¬ 
mains above the observed IGL, the lower limit 
of the EBL, even when DGL is subtracted from 
the sky brightness. As summarized in Hauser & 
Dwek (2001) and Dwek & Krennrich (2013), sev¬ 
eral studies have been modeled the intensity and 
the spectrum of the EBL at redshift z = 0 by dif¬ 
ferent methods (e.g., Stecker et al. 2006, Mazin 
& Raue 2007, Franceschini et al. 2008, Finke 
et al. 2010, Dominguez et al. 2011). Most of 
these results approach the observed IGL (Madau 
& Pozzetti 2000, Totani et al. 2001, Fazio et al. 
2004) and are several times lower than the present 
EBL results. 

To explain the excess diffuse emission reported 
by Matsumoto et al. (2005), the contribution of 
primordial (Pop-Ill) stars were suggested by Sal- 
vaterra & Ferrara (2003). However, Dwek et al. 
(2005b) emphasized that such a large excess is not 
an extragalactic origin since it would have pro- 
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duced a physically unrealistic intrinsic 7 -ray spec¬ 
trum of the blazar PKS 2155-304. Using theoret¬ 
ical constraints on the formation rate of Pop-Ill 
stars, Dwek et al. (2005a) concluded that Pop-Ill 
stars can contribute only a fraction of the EBL 
intensity. This is consistent with the theoreti¬ 
cal contribution of light from Pop-Ill stars, i.e., 
< 0.1nWm~ 2 sr _1 in the near-IR (e.g., Cooray 
et al. 2012a, Inoue et al. 2013, Fernandez & 
Zaroubi 2013). In addition, current EBL con¬ 
straints derived from the 7 -ray observations, as¬ 
suming the different intrinsic spectra of the sources 
(e.g., Dwek & Krennrich 2005, Schroedter 2005, 
Aharonian et al. 2006, Mazin & Raue 2007, Orr et 
al. 2011, Meyer et al. 2012), require the low EBL 
intensity, close to the observed IGL level. Except 
that Guy et al. (2000) allowed the higher upper 
limit of ~ 60nWm” 2 sr _1 at ~ 1 /zm, most of the 
7 -ray constraints on the EBL are inconsistent with 
the present results at 1.25 and 2.2 gm. 

In addition to Pop-Ill stars, several studies re¬ 
cently calculated the ’’exotic” sources’ contribu¬ 
tion to the EBL, such as intrahalo light (IHL), ac¬ 
creting direct collapse black holes (DCBH), and 
dark stars (DS). The IHL could be created by 
tidally stripped stars from their parent galaxies 
by mergers and collisions (Cooray et al. 2012b). 
The IHL intensity estimated by Zemcov et al. 
(2014) is ~ 7 and ~ 2nWm _ 2 sr _1 at 1.25 and 
2.2 /im, respectively. Therefore, the IHL plus the 
observed IGL intensity approaches the EBL re¬ 
sults derived by Wright (1998)-based ZL model 
(Gorjian et al. 2000, Wright 2001, Levenson et 
al. 2007). Cooray et al. (2012b) and Zemcov et 
al. (2014) also suggested that the IHL can ex¬ 
plain the excess in the power spectrum of the dif¬ 
fuse near-IR background, reported by Cooray et 
al. (2012a), Kashlinsky et al. (2005), Cooray et 
al. (2007), Thompson et al. (2007), Matsumoto 
et al. (2011), and Kashlinsky et al. (2012). To ex¬ 
plain the excess in the power spectrum, Yue et al. 
(2013) suggested another candidate, i.e., DCBHs 
in the early universe. The contribution of DCBHs 
to the EBL intensity has a peak at ~ 2 /tm and 
is less than ~ lnWm _ 2 sr _1 at IR wavelengths 
(Yue et al. 2013). Dark stars are the hypothetical 
objects powered by annihilation of either accreted 
or captured weakly interacting massive particles 
before the standard nuclear fusion. Maurer et 
al. ( 2012 ) separately estimated the contribution 


of the colder DS and the hotter ones. As a re¬ 
sult, the contribution of the hotter DS, marginally 
consistent with the current EBL observation, has 
a peak at ~ 2 /im and the intensity is ~ 10 and 
~ 20nWm _ 2 sr _1 at 1.25 and 2.2 /nn, respectively. 
The sum of these exotic sources’ contribution can 
reach the intensity of the present result at 2.2 /trn. 
In contrast, the total of these objects contributes 
less than ~ 20nWm~ 2 sr _1 at 1.25/tm, indicating 
that the present result is approximately two times 
higher than the estimated contribution of the ex¬ 
otic sources plus the observed IGL. 

In conclusion, considering the 7 -ray constraints 
and the currently suggested extragalactic sources’ 
contribution, it is increasingly difficult to attribute 
all of the derived isotropic emission (called “EBL” 
in this paper) to an extragalactic origin, especially 
at 1.25 /tm. Therefore, the excess isotropic com¬ 
ponents may contain light from the local universe; 
the Milky Way and/or the solar system. To iden¬ 
tify where the excess light comes from, we need 
more detailed investigation on the local universe 
as well as extragalactic studies. 

6. SUMMARY 

We reanalyzed the COBE /DIRBE data at 1.25 
and 2.2 /tm. In particular, we measured the EBL 
and evaluated the DGL in the near-IR using the 
DIRBE data, which have wide sky coverage con¬ 
taining regions of both low and high interstellar 
100 /tm intensity. To measure the contribution of 
the starlight at each point in the sky, the ISL 
intensity in each region was calculated from the 
2MASS PSC, which covers almost the entire sky. 

Applying a minimum y 2 analysis at high Galac¬ 
tic latitudes (| 6 | > 35°), we decomposed the sky 
brightness observed by the DIRBE into its four 
components: ZL, DGL, ISL, and isotropic emis¬ 
sion. The DGL was positively linearly correlated 
with the 100 /im brightness, confirming that the 
DGL exists at high Galactic latitudes in the 1.25 
and 2.2 /tm bands. The DGL, which is 1-2 orders 
of magnitude fainter than the other components, 
was extractable because of the high-quality wide- 
held data of the DIRBE. 

The residuals determined by the fitting in¬ 
creased toward the low Galactic latitude region. 
We suggested two possible causes of this phe¬ 
nomenon: faint stars that are filtered out from 
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Fig. 7.— Current measurements of the EBL in the optical and near-IR. Filled red and black circles are 
the results of the present study and Pioneer 10/11 (Matsuoka et al. 2011), respectively. Other EBL 
measurements were obtained from 776T/WFPC2 (Bernstein 2007) (filled black diamonds), IRTS (Matsumoto 
et al. 2005) (filled squares), and AKARI (Tsumura et al. 2013a) (filled triangles). Colored symbols are the 
results of previous C'ORE’/DIRBE data analyzed by Hauser et al. (1998) (pink down-arrows), Cambresy 
et al. (2001) (purple diamonds), Levenson et al. (2007) (green diamonds), Gorjian et al. (2000) (blue 
diamonds), Wright (2001) (cyan diamonds), and Levenson & Wright (2008) (an orange diamond). A brown 
diamond at 3.5 /mi is the result estimated by Dwek & Arendt (1998), assuming the present results at 2.2 fin i. 
Open diamonds, triangles, and squares are the integrated brightness of galaxies obtained from Subaru Deep 
Field (Totani et al. 2001), Hubble Deep Field (Madau & Pozzetti 2000), and Spitzer/IRAC (Fazio et al. 
2004), respectively. For clarity, some results are shifted a little from their exact wavelengths. 
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the 2MASS PSC by the nearest bright stars might 
remain at lower Galactic latitudes or the spatial 
distribution of the ISL might differ between faint 
and bright stars. 

Previous studies have investigated the low 
100 pm region, where the DGL contribution was 
found to be low in the present analysis. After sub¬ 
tracting the ZL, DGL, and ISL, the intensity of 
isotropic emission in our study is approximately 
equal to that of the previous studies. In addition, 
the deviations from isotropy were found to be less 
than 10% in the entire sky at high Galactic lat¬ 
itudes (| 6 | > 35°) at 1.25 and 2.2pm. Although 
the EBL intensity depends on choice of ZL models, 
it shows excess against the observed and expected 
IGL and at both investigated wavelengths. 

Specifically at 1.25 pm, the derived isotropic 
emission is approximately two times higher than 
the observed IGL plus the sum of the contribution 
of suggested extragalactic objects (i.e., Pop-Ill 
stars, intrahalo light, direct collapse black holes, 
and dark stars). In addition, the derived isotropic 
emissions are larger than most of the 7 -ray upper 
limits at both 1.25 and 2.2 pm. Therefore, it is 
possible that the excess emission originates from 
the local universe; Milky Way and/or the solar 
system. 

The ISL maps created from the 2MASS PSC at 
| 6 | > 20 ° are available in the online version of this 
journal. 
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